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The study of excited hadron spectra using Lattice QCD is currently evolving. An important step 
toward obtaining resonance parameters from Lattice QCD is the calculation of finite volume en- 
ergy spectra. Somewhat more rigorous studies of finite volume spectra are currently possible and 
should be completed in the near future. The inclusion of disconnected diagrams is increasingly 
commonplace and the simplest systems which involve mixing between single- and multi-hadron 
interpolating fields are being studied. Advances in all-to-all algorithms which have been crucial 
in this progress are reviewed and a survey of current results is given. Nevertheless, such results 
are preliminary and a thorough discussion of systematic errors is required. We discuss several 
such sources of error, focussing on excited state contamination and the use of the generalized 
eigenvalue problem. Also, the calculation of matrix elements between finite volume Hamiltonian 
eigenstates is discussed. 
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Excited hadrons are interesting to study on the lattice, as their properties cannot be calculated 
using perturbation theory and exhibit many features of strongly coupled QCD. Furthermore, there 
are several unresolved questions related to excited hadrons. For example, most quark model calcu- 
lations of baryon spectra predict considerably more states than have been observed experimentally. 
This may be either because the quark models possess unphysical extra degrees of freedom or be- 
cause most experimental analyses have focused on states which decay to Nn. A summary of the 
current status of this 'missing resonance' problem can be found in Ref. |[l]]. 

In addition to missing resonances, it remains a puzzle why the quantum numbers of known 
hadrons can be predicted from the constituent quark model, which assumes that mesons contain two 
and baryons three quarks. States which have quantum numbers which do not fit these predictions 
are said to be 'exotic' and are the subject of experiments at JLab, Cleo-c, BES III, and PANDA 
(For an overview see Ref. [Q]). 

Lattice QCD calculations are necessarily performed at finite lattice spacing (a) and finite sys- 
tem size (L). Extracting infinite volume resonance parameters from this finite volume data is dif- 
ficult. The spectrum of the QCD Hamiltonian in infinite volume consists of stable single hadron 
states and a continuum of multi-hadron states above threshold. The situation in finite volume is 
quite different; resonance states which were unstable in infinite volume are now eigenstates of the 
finite volume Hamiltonian. Near threshold the finite volume energy eigenvalues are distorted from 
their non-interacting behavior as the spatial extent of the system is varied and an 'avoided level 
crossing' occurs. An analogous situation occurs in the study of string breaking in QCD [||], where 
an avoided level crossing occurs as the length of the string is varied. 

Somewhat more formally, resonances are observed by experiment as rapid variations in the 
cross section as the center-of-mass energy moves through the resonance mass. This can be inter- 
preted as a singularity in the S'-matrix, which (by the LSZ reduction formula) corresponds to a 
singularity in the n-point correlation functions of the underlying quantum field theory. The singu- 
larity structure of a two-point correlation function with suitable quantum numbers is made apparent 
by employing the Kallen-Lehmann representation: 



where A(/7) is the two-point function (the spinless case is shown here) in momentum space, and 
p(/x^) is the spectral density. In infinite volume p{ix^) consists of 5-functions corresponding to 
stable particles below threshold and a continuum of states above threshold. 

Conversely, in finite volume p(/i^) consists entirely of 5-functions corresponding to stable 
states. However, infinite volume scattering phase shifts may be extracted from finite volume lattice 
data below inelastic thresholds. This method |Q] relates infinite volume elastic scattering phase 
shifts to the distortions of finite volume energy spectrum near thresholds mentioned above. 

While elastic scattering phase shifts may be extracted in this manner, systematically extracting 
infinite volume resonance parameters from finite volume data is still an active area of research ||5|, ^, 
^. Regardless, an important first step is to extract the spectrum of QCD in a finite box. This is hard 
enough, as there are many sources of systematic error which much be controlled. As resonances are 
expected to be spatially extended objects, large volumes are needed. Furthermore, extrapolations to 
the physical quark masses are complicated by the inapplicability of chiral effective theory at scales 
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Figure 1: Preliminary results for the chiral behavior of some excited meson states Jl Results for the pion 
channel are shown on the left, while results for the kaon channel are shown on the right. These spectra are 
calculated on ensembles with = at = 0.13fm and Lj ~ 2.2fm, where Aj and a, are the spatial and temporal 
lattice spacings (respectively) and Lj is the spatial lattice extent. Multi-hadron operators are not included in 
analysis, which may result in systematic errors due to threshold effects. 



around and above trip . Additionally, higher energy resonances may suffer from large discretization 
effects so the continuum limit must be taken. Finally, in order to extract the energy of an excited 
state, a large basis of interpolating fields must be employed, and contamination from unwanted 
excited states is present. 

Progress has been made recently toward determining the pion mass dependence of some low 
lying hadron resonances. Current Lattice QCD simulations are performed at unphysically large 
light quark masses and physical results must be obtained by extrapolating to the correct values. 
As discussed above, the use of chiral effective theory to predict the dependence of excited state 
energies on the quark mass is limited. Indeed, if thresholds open up as the pion mass is decreased 
toward the physical point, the chiral behavior of excited states may not even be analytic. Therefore, 
this dependence is largely unexplored and must be determined empirically. Preliminary calcula- 
tions of the chiral behavior of some excited meson states are shown in Fig. [lj, while examples for 
baryons are shown in Fig. |2[ Additionally, a preliminary scan of excited baryon spectra at several 
pion masses can be found in Refs. M, while preliminary results for D mesons and charmionium 



states (without disconnected diagrams) can be found in Ref. [IC]. However, all of these examples 
do not include multi-hadron operators in their analyses, which may introduce threshold-related 
systematic errors. 

We now discuss in more detail the calculation of excited state energies in a finite volume. 
Energies in lattice field theory are typically extracted from the exponential fall-off of temporal 
correlation functions C{t) = (^(f)^(O)) between interpolating fields with particular quantum 
numbers. The energy of the lightest state with these quantum numbers is obtained from the asymp- 
totic (in euclidean time t) single-exponential behavior of the correlation function. Excited states, 
however, can be extracted by forming a correlation matrix Cij{t) = (^,(f)^y(0)) between a set of 
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Figure 2: Preliminary results for the chiral behaviour of some excited baryon states. Results for the odd 
parity nucleon channel are shown on the left [ p^ while states in the A channel [ pj| ] are shown on the right. 
Results were obtained on ensembles with as — at — 0.09fm and Ls — 2.9fm. On the left plot, multi-hadron 
thresholds are indicated. For the A channel, a partially quenched strange quark was used to reproduce the 
physical kaon mass. No multi-hadron operators were included in the analysis, which may result in systematic 
errors due to threshold effects. 



interpolators {^i} and solving the generalized eigenvalue problem (GEVP) [14, 0, 15 1 

C{t)Vn{t,to) = k,r{t,to)C{to)Vn{t,to) (2) 

limE'/^{t,tQ) = lim-dtlnXn{t,to) =E„, 

where is the energy of the nth state interpolated by the {^i}. 

The asymptotic corrections to En^^{t,to) can take different forms based on the relation of tQ to 
t. Generically, it has been proven [Q] that asymptotically En^^{t,to) = £"„ + ^(e^(^"+'^^"^'^), where 
is the energy of the state above En- This correction may be large in systems with closely 



spaced energy levels. However, if the condition to > t/2 is maintained, it has been proven [15] 
that £',j^^(f,fo) =E„ + ff{Q-iEN+i-E„)f^^ where is the dimension of the GEVP. In this manner, the 
corrections can be systematically improved by increasing the size of the GEVP basis. However, this 
may increase the condition number of the resultant correlation matrix, resulting in larger statistical 
errors on the GEVP eigenpairs. 

Due to these two considerations, selecting a GEVP basis is a delicate procedure. As mentioned 
above, increasing the number of operators in the basis may decrease the asymptotic corrections 
but care must be taken to prevent a large condition number A simple procedure to construct a 
basis of operators amounts to applying different levels of Gaussian smearing. However, in the 



case of pseudoscalar static-light mesons, it has been suggested [ ]16[ ] that inteipolating operators 
constructed from different smearings may have poor overlap with excited states. While ideal for 
extracting ground state properties, other types of interpolators may be more suitable for excited 



states. Spatially extended operators designed to transform irreducibly under lattice symmetries [ 17 



T8| , |I9|] have been effective in extracting higher excited states. 

The GEVP may also be used to define 'optimized' interpolating fields [ p^ ] = {vny{t,to)ffi, 
where t and to are typically fixed. Correlation functions of these optimized fields are the diagonal 
elements of the rotated correlation matrix and have increased overlap with the nth state. However, 
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Figure 3: Results from a 3 x 3 GEVP in a toy model consisting of 20 evenly spaced energy levels. Top row: 
The effective energies for the first three levels when the matrix of overlaps is given canonical values. 
Bottom row: The effective energies when the overlaps of all operators with the first level / = 1,2,3) 
are decreased by an order of magnitude. In the case where a low-lying state has small overlaps with all 
operators in the basis, the asymptotic behaviour does not set in until large times. In realistic computations 
statistical errors grow exponentially with time and such non-asymptotic behaviour could be mistaken as a 
'false plateau' . 



mixing can occur and the off-diagonal elements of the rotated correlation matrix must be small to 
extract excited state energies. 

In practice, the maximum temporal separation at which the correlation functions can be evalu- 



ated is limited by the signal-to-noise problem (see e.g. Refs. [ |20| , |21p ) and the finite temporal extent 
of the lattice. It is therefore crucial that Ef/^ {t,to) behaves asymptotically as quickly as possible. 
Although the asymptotic corrections to effective energies in the GEVP are independent of which 
interpolators are included in the set, the time at which asymptotic behavior sets in can vary. In 
particular, if there is a state below £"„ with which the operators have little overlap, the asymptotic 
behavior may set in only at large times. 

To illustrate this point we consider a toy model, which is solved numerically (variants of which 
have been considered in Ref. [p^]). The model is specified by providing an analytic expression for 
the correlation matrix 

C^f {t) =Y^Win.W]m^-^'"' . (3) 

m 

where ro^,,, = m, m = 1...20 and the 3 x 20 matrix y is chosen empirically from an approximate 
calculation of the overlaps in the psuedoscalar static-light meson system [[T^]. The results from the 
3x3 GEVP are shown in Fig. ^ The effective energies corresponding to 'normal' values of the 
Y matrix plateau rather quickly and have corrections of a standard form. When / = 1, ..3 are 
decreased by two orders of magnitude the situation is quite different, however, and the asymptotic 
behavior sets in only at large times. While statistical errors are not included in this model, it seems 
plausible that in the second case with statistical errors one may mistake non-asymptotic behavior 
as a 'false plateau', illustrating the danger of states with small overlaps with all the operators in the 
basis. 
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Figure 4: Types of valence quark line diagrams which must be included for a comprehensive scan of the 
finite volume excited hadron spectrum. Left: Valence quark line connected and disconnected diagrams con- 
tributing to a nucleon-pion correlation function. Right: Connected and disconnected diagrams contributing 
to a flavor singlet meson correlator. 

Given these considerations, a reliable calculation of the hadron spectrum above threshold 
should include both single and multi-hadron interpolating operators (and perhaps other 'exotic' op- 



erators as well [22]). These multi-hadron operators require all-to-all propagators, which describe 
quark propagation from all initial sites to all final sites. This is to be contrasted with 'point-to-all' 
propagators, which describe quark propagation from a single initial site to all final sites. All-to-all 
propagators are also required in flavor-singlet single hadron states and are therefore necessary for 
a comprehensive scan of the meson spectrum. 

The situation can be illustrated by examining valence quark line diagrams. Four such diagrams 
are shown in Fig. ^ When evaluating e.g. a nucleon-pion correlation function both valence quark 
line connected (where all valence quark lines propagate between the initial and final times) and 
valence quark line disconnected diagrams will contribute. Similarly, for flavor singlet mesons both 
connected and disconnected diagrams will contribute. 

Quark propagators are obtained by inverting the lattice Dirac matrix. This matrix is of large 
( 1 2 X {Ls /as )^ X {Lt/ot)) dimension and its inverse cannot be calculated directly. Instead, the action 
of the inverse on a 'source' vector can be obtained by the use of suitable algorithms. Clearly, it is 
not feasible to evaluate all-to-all propagators naively, i.e. by inverting on (12 x [Ls/qsY x [Lf/at)) 
sources, each with support on a single space-time point, spin, and color. 



A novel alternative to naive all-to-all has recently been proposed [|23|]. This 'distillation' 
method calculates quark propagation from a subspace spanned by the low-lying eigenmodes of 
the gauge-covariant Laplace operator to all sites. The projection onto this subspace may be viewed 
as a 'smearing' procedure which creates interpolators with enhanced overlap with the low-lying 
states of interest. Indeed this projection operation preserves all the symmetries of the unsmeared 
propagator and the width of the smearing operator may be controlled by the number of low-lying 
modes contained in the subspace (see Fig. |5]). 

Distillation allows one to compute exact all-to-all propagation from the subspace spanned by 
the low-lying modes of the gauge-covariant Laplace operator and requires a number of Dirac matrix 
inversions proportional to the dimension of this subspace. Unfortunately, the number of modes 
required to maintain a constant smearing radius increases with the spatial volume. As the cost 
of each Dirac matrix inversion is also proportional to the volume, the total cost of this algorithm 
scales like ~ L^. This can be seen by examining the density of the low-lying eigenmodes for several 
volumes, shown in Fig. ^. 

Nevertheless, distillation has been useful in small spatial volumes (L, < 2.5fm). Results from 



6 



Exited hadrons in Lattice QCD 



John Bulava 



0.9 - 
0.8 - 
0.7 - 
0.6 - 
' 0.5 - 
0.4 - 
0.3 - 
0.2 - 
0.1 - 



4 5 

r/a. 



■ 


N=64 


• 


N=32 


• 


N=8 



0.6 
0.5 
.^=0.4 
0.3 
0.2 
0.1, 



T ' 1 ' 1 ' 1 



12 lattice 



,1 „.■" 16^ lattice 

1 I 1 I 1 I 1 L_ 







10 20 30 40 50 
n 



Figure 5: Left: Profile of the smearing operator as a function of the number of eigenvectors on an n/^ = 2 + 1 
lattice with = 3.5flr = 0.12fm and Lj = 2.4fm (see Ref. [^3|]). Right: Low-lying eigenvalue spectra of 
the gauge covariant Laplace operator for two volumes, taken from Ref. [pi]]. Here a, — 3.5flr = 0. 12fm and 
Lj = 1 .4, 1 .9fm. The number of modes contained between the dashed lines increases linearly with the spatial 
volume. 
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Figure 6: Results of a preliminary calculation of the isoscalar meson spectrum (Taken from Ref. 1 28 1) on an 
Hf — 2 + 1 lattice with = 3. So? — 0.12fm, nij^ — 400MeV, and = 1.9fm. Left: Valence quark connected 
and disconnected contributions to an isoscalar scalar correlation function. Disconnected contributions are 
labeled 'D\ while connected contributions are labeled 'C. Right: the effective energy for the lowest states 
in the rj and Tj' meson channels. A 5f = 3 approximation of the temporal derivative is used. 



a preliminary calculation of isoscalar meson spectra using distillation are shown in Fig. ^, while 
preliminary charmonium results which use distillation are shown in Fig. ^ Distillation has also 
been used in calculations of tt — tt scattering phase shifts (see Fig. 1 1 and Refs. [25, 26]) as well as 
preliminary excited baryon and iso vector meson [27] spectrum calculations. 

Clearly, in order to move to large volumes a new all-to-all algorithm must be devised. One 
promising candidate is the 'stochastic LapH'(Laplacian Heaviside) approach [24]. This introduces 
noise in the subspace spanned by the low-lying modes only. Of course this noise may be 'di- 
luted' [3C] in the spin, time, and eigenvector indices. The superiority of these types of diluted 
stochastic sources over conventional dilution is illustrated in Fig. |} 

Furthermore, the volume dependence of the stochastic LapH method is relatively mild for a 
moderate amount of dilution. Formally, this means that the number of dilution projectors (and thus 
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Figure 7: Preliminary results for charmonium spectra (taken from Ref. [^) on an n/ = 2 + 1 lattice with 
Gs = 3.5fl, — 0.12fm, niji — 400MeV, and Ls — 1.9fm. Multi-hadron operators were not included in the 
analysis and non-interacting multi-hadron energies are denoted by dotted lines. 
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Figure 8: Properties of 'stochastic LapH' noise, taken from Ref. [g4|]. Left: The ratio of the error using 
stochastic LapH estimates to the error using the exact (distillation) method for a typical observable (trian- 

— 1/2 

gles), plotted against A^^ ' , where No is the number of dilution projectors. Conventional dilution (squares) 
is also shown for comparison. Results were obtained on a = 2 + 1 ensemble with = 3.5fl, = 0.12fm, 
m^: = 400Me V, and Ly = 1 .9fm. Right: the same quantity for two volumes using stochastic LapH noise and 
dilution. Results were obtained on a ny^ = 2 + 1 ensemble with Oj = 3.5fl, = 0.12fm, m„ = 400MeV, and 
Ls — L9fm( triangles) and 2.4fm (squares). 



the number of inversions) may be held constant as the volume is increased, without degrading the 
quality of the stochastic estimate. The cost now scales as ~ (due to the cost of each inversion) 
rather than the ~ required for exact distillation. This mild volume dependence is demonstrated 
in Fig. I 

Of course, in our estimate of the cost scaling, we have considered the cost of the Dirac matrix 
inversions only. While at moderate lattice volumes (L^ ~ 3fm) this dominates the computational 
cost, the calculation of the Laplacian eigenpairs must be taken into account. This is typically 
done using a variant of the Lanczos algorithm, with some form of polynomial preconditioning. In 
order to ensure numerical stability, a global reorthogonalization of the Lanczos vectors must be 
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Figure 9: Comparison of stochastic LapH and distillation for disconnected diagrams (taken from Ref. [|3l|]) 
on an «^ = 2 + 1 lattice with = 2>.5a, = 0.12fm, m^j =400MeV, and Ls = 1.9fm. Results using distillation 
are shown with black circles, while the ideal stochastic LapH dilution scheme is shown with with blue 
triangles. The left plot shows the 'box' diagram contribution to / = n — n scattering while the right 
plot shows the disconnected contribution to a / = scalar correlation function. In the case of the scalar, 
the stochastic LapH dilution scheme requires 1024 inversions per configuration while the distillation result 
requires 16384. 



performed periodically. This reorthogonalization alone scales like so at larger volumes the cost 
of generating the Laplacian eigenpairs may become significant. 

Apart from the valence quark line connected observable shown in Fig. ^ this method is also 
adequate for the estimation of disconnected diagrams. Flavor singlet correlation functions con- 
taining disconnected diagrams have a severe signal-to-noise problem, as the variance contains a 
component which is independent of the temporal separation. Therefore, even the variance of the 
exact all-to-all result is large, which allows one to reasonably estimate these diagrams using a mod- 
erate amount of dilution. Indeed, in Fig. ^ it is shown that an error similar to the exact distillation 
result can be obtained stochastically with e.g. a factor ~ 16 fewer Dirac matrix inversions per 
configuration. 

After tests confirmed its utility, the stochastic LapH technique was applied to a preliminary 
calculation of the meson spectrum in a larger volume, the results of which are shown in Fig. [I^. 
Although multi-hadron operators are not yet included, such a calculation would not have been 
feasible without the stochastic LapH approach. 

Of course, there has been work on the inclusion of multi-hadron operators as well. The sim- 
plest system which exhibits mixing between single- and multi-hadron interpolating fields is the 
7 = 1 K — K ip) sector. Several studies of the p -meson [33, 34, 36] have been presented at 
this conference. A summary of the unquenched results as well as results from a smaller volume 



calculation of the 7 = 1 tt — tt phase shift using distillation are shown in Fig. 1 1 



Flavor singlet single hadron correlation functions also require valence quark line disconnected 
diagrams and can similarly be computed. Not only are flavor singlet meson spectra interesting 
in their own right, but may also appear in other channels as multi-hadron decay product states. 
To this end, flavor singlet interpolating operators must be studied both at rest and with non-zero 
spatial momentum if they are to be included as multi-hadron states in a GEVP analysis. Some 
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Figure 10: Isovector meson results (taken from Ref. [^) on a n/^ = 2 + 1 ensemble with cis = 3.5a, = 
0.1 2fm, niji = 400MeV, and Lj = 2.9fm. The horizontal axis denotes different irreducible representations 
(uTeps) of the lattice symmetry group. The subscript 'g' denotes even parity irreps, while 'u' denotes odd 
parity. The '+' and '-' labels refer to the G-parity of the irrep. 7 = states are found in the Aj irrep, J — I 
in the Ji, 7 = 2 in the E and T2 irreps, and 7 = 3 in the A2, T\, and T2 irreps. Multi-hadron operators are not 
included in the GEVP analysis and multi-hadron thresholds for each set of quantum numbers are indicated 
by the shaded area. 
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Figure 11: Results on the p meson from unquenched ensembles. Left: A compilation of recent unquenched 
results taken from Ref. | ]33[ | (color online). The black points are taken from Ref. [^, the red points from 
Ref. and the open circle from Ref. Riglit: The I^\n — K scattering phase shift using distillation 
(taken from Ref. on an = 2 ensemble with Oj = a, = 0.12fm, m,i ~ 260MeV, and Ly = 1 .9fm. 
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Figure 12: Preliminary results which require valence quark line disconnected diagrams using the stochastic 
LapH method. Top Row: Effective energies from a 5x5 GEVP in the Isoscalar pseudoscalar sector (rj 
mesons) with one unit of lattice momentum, taken from Ref . 1 38 1 . Results are from an rif = 2+1 ensemble 
with flj = 3.5fl, = 0.12fm, m^i = 400MeV, and Lj = 1.9fm. Bottom Row: Effective energies from a 4x4 
GEVP in the isoscalar scalar sector (taken from Ref. [|3^]) on the same ensemble. The basis of operators 
used in this analysis consists of a local a-meson operator, two I = n — n meson operators, and a scalar 
glueball operator constructed from the low-lying eigenvalues of the gauge-covariant Laplace operator. The 
mixing between all three types of operators is significant. 



preliminary results for moving rj mesons are shown in Fig. |l^. Such a calculation is complicated 
by the reduced symmetry group of particles in motion and operators which transform according to 
irreps of the lattice little group must be constructed for various lattice momenta [37]. 

Fig. |l^ also shows a first glimpse at a more realistic spectrum calculation in finite volume, 
specifically the isoscalar-scalar (vacuum) sector. This channel is particularly difficult as glueball 
states are present, in addition to single and multi-hadron states. A realistic calculation of the spec- 
trum in this channel therefore requires interpolating fields with a reasonable overlap to these three 
types of states. Results of a preliminary such calculation are shown in Fig. |l^. There a GEVP 
analysis was performed which contained two local (a-meson) operators, two I = n — n operators 
with equal and opposite momenta, and a scalar glueball operator created using the low-lying eigen- 
values of the gauge covariant Laplace operator. Indeed, significant mixing was found between all 
three types of operators in this channel. 

Apart from finite volume energies, matrix elements of local operators between finite volume 
Hamiltonian eigenstates can also be calculated. In general, it is non-trivial to rigorously relate 
these quantities to those which have a well-defined infinite volume limit, but nonetheless these 
matrix elements may have phenomenological implications. Apart from the work presented here, 
preliminary calculations of transitions between excited and ground states (which neglect the effect 
of disconnected diagrams) have been performed in the nucleon and charmonium sectors [39, 4C]. 

The situation is somewhat simplified for matrix elements of the light-light axial current be- 
tween B and B* mesons. First, no disconnected diagrams contribute to the required three-point 
correlation functions. Secondly, radial excitations of pseudoscalar static light mesons can only de- 
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Figure 13: Effective matrix elements of the Ught-Hght axial current between radial excitations of B* and 



B mesons in the static limit. These results (taken from Ref. |43|) are from an nj — 2 ensemble with 
Oj = a, = 0.065fm, nij^ = 447MeV, and — 2.1fm. Although the renormaUzation factors for this lattice 
discretization are known non-perturbatively, only bare matrix elements are shown here for illustrative pur- 
poses. Nonetheless, since the axial current is taken here at zero momentum transfer, these matrix elements 
are multiplicatively renormalized only and ratios between bare matrix elements are renormalized quantities. 



cay strongly via the emission of two pions. At pion masses which are large enough so that the first 
radial excitation is significantly below the two-pion threshold, this excitation is stable and matrix 
elements of the axial current involving this state and the ground state have a well defined infinite- 



volume limit. Results from a preliminary calculation of these matrix elements are shown in Fig. 13 



where rather than effective energies, effective matrix elements are plotted. These effective matrix 



elements [[151, |16p are also obtained from solutions of the GEVP and are defined as 

M'£\t) = R,„{t,to)Rn{t,to) X {v,„{t,to),C^'"{t,to)v„{t,to)), (4) 

where parentheses denote an inner product over the GEVP indicies, v„,(f,fo) are GEVP eigenvec- 
tors, and Rfn{t,to) are normalization factors constructed to cancel the asymptotic time dependence. 
Unlike the effective energies, here the condition to > t /2 is not required to reduce the asymptotic 
corrections, but t must be at least larger than tQ. The asymptotic behavior of this effective matrix 
element is proven [Oil to beM^J (f,fo) = {B* ,m\Ai\B,n) + 0'(e^iEN+i-E„,.„)to-^^ where A,- is the light- 



light axial current, while and \B,n) are finite volume Hamiltonian eigenstates corresponding 

to radial excitations of static-light zero-momentum B* and B mesons, respectively. 

The required matrix of three point correlation functions is given as C^j" {t) = {i)'^{2t)Ak{t)S'{0)) 
where and {i^,} are sets of interpolating operators for B* and B mesons, respectively. Since 
the three-point correlation function contains two separations which both must be taken large, the 
condition to < t means that in terms of the total separation = 2t the asymptotic correction to 
M^J{t,tQ) is i^(e-(^'v+i-£«,„)fi/2)^ compared with the ff(^e^iEN+i-E„)t,^^ discussed earlier for effec- 
tive energies. 

It should be noted that for matrix elements for which E,„ = E„, the asymptotic corrections may 
be improved [ffTl W^, 42|. This entails the use of a summed insertion where instead of C^P'{ti,t2), 



we employ D^P'{t) = C^p' {t + t\,ti). Using this summed insertion with the GEVP it is possible 



to obtain an asymptotic correction which is proven [16] to be i^(e^(^'^+'^^'')''), as in the effective 
energies. However, effective matrix elements constructed from this summed insertion are typically 
noisier, as a temporal derivate must be taken numerically. 

In conclusion, excited hadron spectroscopy in Lattice QCD is a field which is currently evolv- 
ing. While the systematic extraction of infinite volume resonance parameters from finite volume 
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energy spectra remains a difficult problem, progress has been made in the extraction of finite vol- 
ume energy spectra from lattice data. One of the major sources of difficultly in calculating finite 
volume energy spectra in Lattice QCD is the need to include both single and multi-hadron inter- 
polating fields in a correlation matrix analysis. The requires the evaluation of valence quark line 
disconnected diagrams which require the knowledge of quark from all initial space-time points to 
all final space-time points. 

Although all-to-all propagators cannot be calculated naively, an efficient stochastic algorithm 
(stochastic LapH) has been developed which requires a computational cost that scales linearly 
with the spatial volume. This technique (as well as others) has enabled preliminary calculations 
of spectra in systems which require disconnected diagrams, such as the 7 = scalar sector and 
the p -meson sector, where multi-hadron diagrams should be included in the basis of interpolating 
fields. Finally, apart from finite volume energies, transitions between finite volume Hamiltonian 
eigenstates can be calculated. 
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